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MULTISIGNAL  MINIMUM-CROSS-ENTROPY  SPECTRUM  ANALYSIS 

WITH  WEIGHTED  PRIORS 


INTRODUCTION 

Multisignal  minimum-cross-entropy  spectrum  analysis  (MCESA)  is  a  method  for  estimating  the 
power  spectrum  of  one  or  more  signals  when  a  prior  estimate  for  each  is  available  and  new  information 
is  obtained  in  the  form  of  values  of  the  autocorrelation  function  of  their  sum  [1],  The  resultant  esti¬ 
mates  are  the  solution  of  a  constrained  minimization  problem:  they  are  consistent  with  the  autocorrela¬ 
tion  information  and  otherwise  as  similar  as  possible  to  the  respective  prior  estimates  in  a  precisely 
defined  information-theoretic  sense.  Multisignal  MCESA  [1]  is  a  generalization  of  MCESA  [1,2]  and 
reduces  to  it  in  the  special  case  when  the  number  of  signals  is  one.  Multisignal  MCESA  can  be  derived 
[1,2]  as  an  application  of  the  principle  of  minimum  cross  entropy  [3-5]  or,  alternatively,  by  minimizing 
a  sum  of  Itakura-Saito  distortions  [6]. 

Multisignal  MCESA  applies  when,  for  instance,  one  obtains  autocorrelation  measurements  for  a 
signal  corrupted  by  independent  additive  noise,  and  one  has  some  prior  knowledge  concerning  the  spec¬ 
tra  of  both  the  uncorrupted  signal  and  the  noise.  The  results  are  posterior  signal-  and  noise-spectrum 
estimates  that  take  both  the  prior  estimates  and  the  autocorrelation  information  into  account. 

The  multisignal  MCESA  procedure  presented  in  Ref.  1  gives  the  same  weight  to  each  of  the  prior 
spectrum  estimates— they  are  all  treated  on  the  same  footing.  However,  in  situations  that  arise  in  prac¬ 
tice,  one  may  have  more  reliable  information  about  the  spectra  of  some  of  the  signals  than  about  the 
others.  Consider  a  speech  signal  corrupted  by  additive  background  noise;  suppose  the  background  is 
more  nearly  stationary  than  speech.  If  it  is  possible  to  detect  pauses  in  the  speech  reliably,  then  meas¬ 
urements  of  the  sum  signal  during  a  speech-free  interval  will  yield  an  estimate  of  the  noise  spectrum 
that  can  serve  as  a  prior  noise  estimate  during  an  interval  when  speech  is  present.  The  result  may  well 
be  a  better  prior  estimate  for  the  noise  spectrum  than  any  prior  estimate  that  can  be  obtained  for  the 
speech  spectrum.  Alternatively  one  might  be  able  to  obtain  a  good  estimate  of  the  noise  power  spec¬ 
trum  by  conventional  spectrum  analysis  of  a  signal  from  a  microphone  exposed  to  the  noise  but  not  the 
speech.  In  both  cases  it  would  be  desirable  to  give  greater  weight  to  the  noise  prior  than  to  the  speech 
prior  in  deriving  posterior  spectrum  estimates.  In  other  situations  it  might  be  desirable  to  rely  more 
heavily  on  the  speech  prior  than  on  the  noise  prior. 

It  is  furthermore  possible  to  have  prior  information  about  the  spectrum  of  an  individual  signal 
that  is  more  reliable  in  some  frequency  ranges  than  in  others.  Thus  in  some  situations,  for  example, 
one  might  even  wish  to  give  greater  weight  to  the  noise  prior  at  high  frequencies  and  to  the  speech 
prior  at  low  frequencies. 

We  present  in  this  report  a  generalization  of  multisignal  MCESA  that  allows  a  frequency- 
dependent  weight  to  be  attached  to  each  prior  estimate.  Aside  from  the  weights,  inputs  to  the  pro¬ 
cedure  are  the  same  as  to  multisignal  MCESA:  a  prior  spectrum  estimate  for  each  signal  and  autocorre¬ 
lation  values  for  the  sum.  The  results  again  are  posterior  spectrum  estimates  that  are  consistent  with 
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the  given  autocorrelation  information.  When  all  of  the  weights  are  constant  and  equal,  the  results  are 
identical  to  those  from  multisignal  MCESA.  In  general,  increasing  a  weighting  parameter  tends  to  bring 
the  corresponding  posterior  spectrum  closer  (in  the  sense  of  Itakura-Saito  distortion)  to  the  correspond¬ 
ing  prior  spectrum  at  the  expense  of  increasing  the  distortions  for  other  prior-posterior  pairs.  This  is 
demonstrated  in  the  fourth  section  of  this  report,  where  it  is  also  shown  that  the  posterior  spectra  can 
be  obtained  by  constrained  minimization  of  a  weighted  sum  of  Itakura-Saito  distortions.  In  the 
remainder  of  this  section,  we  describe  the  multisignal  MCESA  method  of  111.  In  the  second  section  of 
this  report,  we  present  the  generalized  method  and  the  heuristic  arguments  that  first  led  us  to  the  gen¬ 
eralization.  In  the  third  section  we  give  a  derivation  based  on  minimizing  a  weighted  sum  of  cross 
entropies,  and  in  the  fifth  section  we  rework  a  numerical  example  from  Ref.  1. 

Multisignal  MCESA  estimates  the  power  spectra  S,(f)  of  a  number  K  of  independent,  real,  band- 
limited,  stationary  processes  (signals)  with  bandwidth  W,  given  a  prior  estimate  P,  for  each  Sh  and 
given  in  addition  the  values  -  R'ol(tr)  of  the  autocorrelation  function  of  the  sum  of  the  processes 
at  finitely  many  lags  t„  r  —  0  ,  ....  M.  The  prior  estimates  Pt  may  be  thought  of  as  the  best  guesses 
at  S,  we  could  make  in  the  absence  of  the  autocorrelation  data.  Under  the  assumption  of  indepen¬ 
dence,  we  write  R,ot  as  a  sum  of  autocorrelation  functions  Rt  of  individual  processes: 

RMU)  -  *,(/), 

where 

R,U)  -  2 fQ  S,(f)  cos  2irfl  df. 

The  multisignal  MCESA  estimator  has  the  form  [1] 

0,00  -  — - *-! -  (1) 

TJJT  +  J) 2Pr  cos  2irftr 

where  the  f}r  are  chosen  so  that  the  Q,  are  consistent  with  the  given  autocorrelation  values: 

Q,(f)  cos  lit  ft,  df.  (2) 

The  parameters  /3,  are  Lagrange  multipliers  that  arise  in  the  solution  of  a  minimization  problem  with 
Eq.  (2)  as  constraints.  The  minimization  problem  can  be  formulated  as  an  application  of  the  principle 
of  minimum  cross  entropy  or,  alternatively,  as  the  minimization  of  the  sum 


of  Itakura-Saito  distortions. 

THE  METHOD 

To  use  the  multisignal  MCESA  estimate  with  weighted  priors,  we  must  supply  not  only  autocorre¬ 
lation  values  R,t0'  for  the  sum  of  the  signals  and  prior  estimates  P,  for  the  individual  signals,  but  also  a 
weight  w,  associated  with  each  P,.  The  w,  may  be  simply  K  positive  constants  or,  more  generally,  K 
functions  w,(f)  of  frequency.  The  estimate  has  the  form 


where  the  parameters  fi,  are  to  be  chosen  so  that  the  constraints  of  Eq.  (2)  are  satisfied. 
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We  first  arrived  at  something  like  Eq.  (4),  with  frequency-independent  weights,  by  considering  a 
somewhat  artificial  situation.  Suppose  that  each  signal  s,  is  the  sum  of  n,  independent  signals  su  with 
spectra  Sy,  and  that  for  each  i  we  have  equal  prior  estimates  PJnt  for  all  the  Su,  j—  l,  ,  n,.  Then 
we  might  consider  nt  to  be  a  measure  of  the  quality  or  reliability  of  P,  as  an  estimate  of  S,  —  £uS0’  by 
analogy  with  the  fact  that  the  sum  of  n  independent  random  variables  distributed  identically  with  x/n 
has  a  smaller  variance  (by  a  factor  of  n)  than  a  single  random  variable  x.  The  MCESA  posterior  esti¬ 
mates  Qy  for  Sy  in  this  situation  are  given  by 


Qij(f)  - 


_ 1 _ 

t>(7)/«;  +  ,5  2f*r  cos  2irA 


where  the  /3,  are  chosen  to  satisfy  constraints 

' QtJ(f)  cos2n ftr  df  (5) 

i-l  J-\ 

We  consequently  obtain 

Q,<J)  -  I  Qu(f)  -  n^if) 

j- 1 


- j - j — H - -  vo/ 

W  +  T,%  2A  cos  2irA 

as  a  reasonable  posterior  estimate  of  S,(f).  When  Q,  is  given  by  Eq.  (6),  the  constraints  of  Eq.  (5) 
become  equivalent  to  those  of  Eq.  (2) . 

The  posterior  spectra  Q,  are  not  altered  if  the  n,  are  replaced  with  a  proportional  family  of 
numbers  wh  since  a  common  factor  can  be  absorbed  into  the  /3,.  Furthermore,  there  does  not  seem  to 
be  much  point  to  requiring  the  ratios  of  the  w,  to  be  rational— we  regard  the  model  of  j/  as  a  sum  of 
components  s0  as  being  suggestive,  but  not  necessarily  to  be  taken  literally.  We  therefore  replace  n,  in 
Eq.  (6)  with  arbitrary  positive  numbers  w,  and  obtain 

«</>  -  -j - ttJ -  m 

7m  *  Til ™  u/l’ 


Allowing  the  weights  to  be  frequency-dependent  can  be  motivated  heuristically  by  considering 
another  somewhat  artificial  procedure— converting  a  single-signal  problem  into  a  multisignal  problem. 
Assume  we  are  given  autocorrelation  values  Rr  and  a  prior  estimate  P  for  the  spectrum  S  of  a  single 
signal  s.  Suppose  we  arbitrarily  partition  the  band  of  frequencies  from  0  to  IV  into  two  bands  B\  and 
B2  and  write  s  as  a  sum 


s(r)  -  sj(r)  +  s2(t) 


of  two  signals  whose  spectra  St  and  S2  are  confined  to  the  respective  bands  B\  and  B2 


S((/) 


s(f),  f  e  Hi 

o.  fKBr 


r  -  1,2.  We  solve  the  two-signal  MCESA  problem  with  autocorrelation  values  given  by  R,tot 
prior  estimates  P2  given  by 


Pi  if) 


Pif),  /*B, 

0,  /  M  ' 
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and  frequency-independent  weights  wf.  Strictly  speaking,  Eq.  (4)  is  ill-defined  where  /*,(/)  —  0;  how¬ 
ever,  there  is  no  problem  with  the  following  alternative  form, 

W) 


(?/(/)  - 


1  +  —P,(f)  £  2/9,  cos  2  IT  ft, 


The  sum  Q  of  the  posteriors  satisfies 

Qif)  -  Qiif)  +  Qi<J > 
|Oi (A  /€ 

“  (&</).  /  € 

and  therefore 


7T7T  +  ^  ,5  2^r  cos  2nftr 

when  f€Bh  (/  —  1,  2).  With  equal  weights  W|  —  wj-  1,  we  thus  merely  recover  the  single-signal 
MCESA  result  in  a  roundabout  way.  However,  we  may  write  this  as 


<200- 


where,  with  general  weights. 


TUT  +  WF  Jf  2Pr  C°S  2nf‘r 


,  „  J»i.  /€  B\ 

w:.  /€*, 


We  could  equally  well  have  partitioned  the  frequencies  into  any  number  of  sets  B„  not  merely  2;  from 
that  point  it  is  an  obvious  step  to  consider  Eq.  (8)  with  weights  w  (/)  not  restricted  to  a  finite  number 
of  values  wt.  We  are  thus  led  to  a  single-signal  version  of  Eq.  (4).  The  generalizations  that  lead  from 
the  unweighted  multisignal  MCESA  estimate  Eq.  (1)  to  Eq.  (7)  and  Eq.  (8)  combine  to  yield  Eq.  (4)  in 
full  generality. 

DERIVATION 

We  collect  here  some  notation  and  results  from  Ref.  1.  For  the  K  signals,  we  use  discrete- 
spectrum  approximations 

N 

S/it)  -  ]£  ia,k  cos  2nfkt  +  b,k  sin  2wfkt), 

k-\ 

(i  -  1  ,  ...  ,  K).  The  /*  are  nonzero  frequencies,  not  necessarily  uniformly  spaced,  and  the  alk  and 
blk  are  random  variables.  We  define  random  variables 


Xu,  -  +  b}k) 


representing  the  power  of  process  s,  at  frequency  fk,  and  we  consider  their  joint  probability  density 

qf(x),  where  x  -  Gq,  .  .  .  ,xk)  and  x,  -  Cx, i . xlN).  Assuming  prior  estimates  Plk  -  Pt(fk)  of  the 

power  spectra  of  the  s„  we  write  prior  estimates  p  of  qf  in  the  form 


p(x)  - 


n  i* 

M 


Ptk(x,k). 
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where 

/  \  1  xik 

Pik'xik '  —  ”J5 —  —p - 

"  ik  "  ik 

Autocorrelation  values  Rrtot  =  R(tr)  for  the  sum  ot  the  s,  are  assumed  known  at  lags  r, 
(r  —  0 . M)  with  /0  =  0.  (In  Ref.  I .  /?r  is  written  in  place  of  Rrroi.)  We  write  these  as  linear  con¬ 

straints 

Rrto1  =  I  I  J CrifcJC^x)  dx  (9) 

/-I  *-l 

on  expectation  values  of  q  \  where 

crk  =  2  cos  2 t rfktr. 

We  obtain  a  posterior  estimate  q  of  qf  by  minimizing  the  cross  entropy 


H(q,  p)  =  f  q  ( x )  log  dx 

subject  to  the  constraints  (Eq.  (9)  with  q  in  place  of  q1)-,  the  result  has  the  form 

K  N 

«(*)-nn  «*(*<*)• 

l-k  k-l 

where  the  qtk  are  related  to  the  posterior  estimates 

Qik  “  Qi  (/*)  -  /  X uA  (x)  dx 

of  the  power  spectra  of  the  s,  by 


Qik(xlk ) 


exp 


~xik 

TZT 


We  find  a  discrete-frequency  version  of  Eq.  (1), 


wnere  the  /3,  must  be  chosen  so  that 


is  satisfied. 


Qik  — 


1 


1  M 

T  +  I  PrC* 

^ik  r- 1 


K  N 

IZ  crkQlk“  Rr°' 

/-I  *-l 


(10) 


(11) 


(12) 


Equation  (10)  states  the  posterior  independence  of  the  xlk.  We  would  have  obtained  the  same 
results  Eq.  (11),  Eq.  (12)  if  we  had  assumed  Eq.  (10)  from  the  start,  choosing  the  qik  to  minimize 
H(q,  p)  subject  to  the  constraints  with  q  expressed  in  the  form  of  Eq.  (10).  For  such  densities  q ,  we 
have 

H (q.  p)  -  jj  Jj  H(qlk,  plk),  (13) 

and  the  constraints  assume  the  form 

I  l/  «■***«*(**)  dxik  “  R'01-  (14) 

/“  *-i 
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Thus  multisignal  MCESA  (without  weighting  factors)  can  be  obtained  by  minimizing  Eq.  (13)  subject 
to  the  constraints  of  Eq.  (14).  Our  generalization  is  to  replace  the  right-hand  side  of  Eq.  (13)  with  a 
weighted  sum 

K  N  K  n  r  aJxJ 

Z  Z  w*  Pik )  “  Z  Z  H'«J  9tk(xik)  log  — r  dxik. 

/-I  *-l  l-l  k-l  Plk'xlk' 

This  is  to  be  minimized  with  respect  to  variations  in  the  q,k  subject  to  the  constraints  Eq.  (14)  — 
together,  of  course,  with  the  normalization  constraints 

fqlk(xlk)dxlk-\.  (15) 

The  result  is 


?/*(*/*)  “  Pik(x,k)  exp 


-  1  - 


^lk 

~*tk 


1  N 
—  Z  PrCrk*,k 
wUc  r-1 


where  the  /3,  and  the  klk  are  Lagrange  multipliers  corresponding  to  the  constraints  Eqs.  (14)  and  (15) 
From  this  it  follows  that  Eq.  (11)  holds  with 

Qik-  f  Xikqlk(xlk)dxlk 

1 

-  - — a - , 

+  J_Z/w 

Wlk  1 


(16) 


the  argument  is  the  same  that  established  (13)  and  (14)  from  Eq.  (11)  in  Ref.  2. 

Passing  to  the  continuous-frequency  case,  we  write  the  weights  as  w,(f)  and  replace  Eq.  (16)  with 
Eq.  (4). 

PROPERTIES 


We  begin  with  two  trivial  observations.  The  first  is  that,  as  noted  in  the  second  section  of  this 
report,  the  weights  wt  may  be  scaled  by  a  common  factor  without  affecting  the  results.  The  second  is 
that  the  present  method  does  indeed  reduce  to  multisignal  MCESA  when  the  w,  are  all  constant  and 
equal. 


Next  we  show  that  the  posterior  spectra  Q ,  can  be  obtained  by  minimizing  the  sum 

of  weighted  Itakura-Saito  distortions  subject  to  the  constraints  of  Eq.  (2).  We  form  the  expression 
k  rw 


(17) 


&(/)  (?/(/)  1  M  K  cw 

~PJ7T  ~  l0g  7,(7)  “  1  \df  +  Xj  2^o  GfC/")  cos  lit  ft,  df 

involving  Lagrange  multipliers  0,,  and  we  set  its  variation  with  respect  to  £),  (/)  equal  to  zero 

I  1  -  1 

W  wi 


»,(/) 


M 


+  £  2/3 ,  cos  lit  ft,  —  0. 


r-l 


This  implies  Eq.  (4).  We  obtain  a  minimum  of  Eq.  (17)  since  the  second  variation,  *v,(/)/0,(/)2,  is 
positive. 


6 
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Finally,  we  justify  the  claim  that  increasing  one  of  the  weights  tends  to  decrease  the  distortion 
between  the  corresponding  prior  and  posterior  spectra  while  increasing  the  distortions  for  other  prior- 
posterior  pairs.  We  write  D(Qt,Pt)  for  the  Itakura-Saito  distortion  of  P,  with  respect  to  Qt  (the  integral 
in  Eq.  (3)).  We  first  consider  frequency-independent  weights.  Let  w\  be  the  result  of  increasing  w,  for 
a  particular  value  a  of  i  and  leaving  the  rest  of  the  weights  the  same: 


H >'a  >  Wa, 

w'j—Wj  (i  ^  a). 

Let  the  use  of  weights  w',  result  in  posterior  spectra  Q\.  We  will  show  that 

D{Q'a,Pa)  <  D{Qa,P„), 

and  that 

DCQ'^P,)  >  D(Q„P,) 

for  at  least  one  value  of  i  (necessarily  different  from  a). 


(18) 

(19) 


Now  Q,  minimizes 

X  w,D(Q,,P,) 

i 

subject  to  the  constraints  (2)  while  Q\  minimizes 

£  w'/DiQ'i.P ,) 

i 

subject  to  the  same  constraints  (with  (?',  in  place  of  Q,).  It  follows  that 

£  w,D  (Q„P,)  <  £  W'DiQ^P,)  (20) 

i  i 

and 

£  w'iD(Qi,P,)  >  £  w'tDtQ'^P).  (21) 

/  t 

Subtracting  Eq.  (20)  from  Eq.  (21)  yields 

(w'a  -  »a)D(Qa,Pa)  >  (w'a  -  wa)D(Q'0,Pa). 

Since  (w'a  -  w„)  is  positive,  we  have  Eq.  (18).  But  it  follows  from  Eq.  (20)  that  Eq.  (19)  holds  for 
some  L 


A  similar  argument  establishes  a  somewhat  similar  result  for  frequency-dependent  weights.  For 
simplicity  we  state  this  for  the  single-signal  case.  Let  »•'  be  the  result  of  increasing  won  some  band  B 
of  frequencies: 

w’(/)  >  "(A  f*B 

*'(/)-"(/),  MB. 


Let  the  use  of  wand  w'  result  in  posterior  spectra  Q  and  Q\  respectively.  Then 


for  some  B'  £  B,  while  the  reverse  inequality  holds  for  some  B"  disjoint  from  B'.  That  is,  increasing 
the  weight  on  B  decreases  the  distortion  on  some  subset  B'  of  B  while  increasing  the  distortion  else¬ 
where. 
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EXAMPLE 


in  this  section  we  take  up  a  two-signai  numerical  example  that  was  used  in  Ref.  1  to  illustrate 
multisignal  MCESA.  We  use  the  same  prior  spectra  and  the  same  autocorrelation  data,  and  we  show 
plots  of  the  posterior  spectra  that  result  from  various  choices  of  the  weights.  The  autocorrelation  data 
were  derived  from  the  sum  of  a  pair  of  assumed  original  spectra  and  SB.  (In  Ref.  1  the  indices 
stood  for  signal  and  background )  These  spectra  each  have  a  sharp  peak  at  a  single  frequency.  They  are 
shown  in  Fig.  1.  The  priors  Ps  and  Ps  were  taken  to  be  respectively  flat  and  identical  to  SB.  They  are 
shown  in  Fig.  2,  which  corresponds  to  Figs.  3  and  2  of  Ref.  1.  Figures  3-7  show  poste  •  spectra  Qs 
and  Qb  for  five  choices  of  constant  weights:  the  ratios  w$fwB  are  90/10,  75725,  50/  -.5/75,  and 
10/90.  The  spectra  in  Fig.  5,  the  equal-weight  case,  are  identical  to  those  in  Figs.  6  and  r  Ref.  1.  In 
each  case  there  is  a  sharp  peak  in  Qa ,  corresponding  to  the  peak  that  is  in  S8  and  PB  each  case 

there  is  a  second  peak  in  one  or  both  of  the  posteriors,  corresponding  to  the  peak  that  Ss  but  in 

neither  prior;  its  appearance  in  the  posteriors  is  due  entirely  to  information  contained  ir  ■'tocorre- 
lation  data.  In  Fig.  3,  for  which  the  prior  Ps  is  heavily  weighted,  the  posterior  Qs  is  app  .%•  ely  flat, 
resembling  the  prior,  and  most  of  the  power  in  the  second  peak  is  attributed  to  Qb  In  F,*,.  ,  PB  is  the 
heavily  weighted  prior;  most  of  the  power  in  the  second  peak  is  attributed  to  Qs,  whose  corresponding 
prior  received  a  small  weight  indicating  prior  uncertainty.  The  figures  between  show  a  progression  from 
the  one  extreme  to  the  other  as  the  weight  shifts. 
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Fig.  I  —  Assumed  original  spectra  for  example 
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Fig.  4  —  Posterior  spectral  estimates  (ws/wB  —  75/25) 
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Fig.  5  —  Posterior  spectral  estimates  (*>§/  wB  -  50/50) 
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Fig.  6  —  Posterior  spectral  estimates  (ws/wB  -  25/75) 
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Fig.  7  —  Posterior  spectral  estimates  (ws/wB  -  10/90) 


DISCUSSION 

The  generalization  of  multisignal  MCESA  we  have  presented  allows  one  to  take  into  account  the 
possibility  that  prior  spectrum  estimates  are  not  equally  reliable  for  all  signals  or  in  all  frequency  ranges 
for  the  same  signal.  One  can  express  greater  or  less  confidence  or  "degree  of  belief  in  the  prior  esti¬ 
mates  for  the  various  signals  in  various  frequency  ranges  by  assigning  greater  or  smaller  values  to  the 
corresponding  weights.  Roughly  speaking,  increasing  a  weight  tends  to  bring  the  corresponding  poste¬ 
rior  closer  to  the  prior.  This  statement  is  made  more  precise  in  the  results  proved  in  the  fourth  section 
of  this  report  and  is  illustrated  in  the  examples. 

Computer  programs  for  multisignal  MCESA  spectrum  estimation  with  weights  may  be  found  in 
the  appendices  of  Ref.  7. 

For  a  single  signal,  somewhat  different  considerations  have  led  Chu  and  Messerschmitt  [8]  to  the 
idea  of  minimizing  a  weighted  Itakura-Saito  distortion.  However,  their  procedure  involves  minimizing 
the  integral  in  Eq.  (17)  by  varying  P  rather  than  Q.  Chu  and  Messerschmitt  view  the  frequency- 
dependent  weight  as  a  means  for  taking  into  account  the  varying  perceptual  importance  of  various  fre¬ 
quency  ranges  of  a  speech  signal.  They  consider  the  weighted  distortion  between  a  true  spectrum  and 
an  all-pole  estimate,  rather  than  that  between  a  posterior  estimate  and  a  prior  estimate,  and  they  thus 
obtain  a  new  method  for  choosing  all-pole  estimates  of  the  usual  form,  rather  than  an  estimate  in  a  new 
form  such  as  Eq.  (4). 
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